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I. INTRODUCTION 


Soliton and black hole solutions of Einstein-Yang-Mills (UYM) theory have been studied 
extensively for over twenty years (see, for example, Ref. Ill for a review). The hrst solu¬ 
tions found were spherically symmetric, purely magnetic, asymptotically flat, solitons^ and 
black holes^ in four-dimensional su(2) EYM. Discrete families of solutions were found nu¬ 
merically and their existence was later proven (see Refs. 4 and 0 for some analytic work). 
The purely magnetic gauge held is described by a single function u, which has at least one 
zero. The families of solutions are characterized by the event horizon radius Vh (with r/j = 0 
corresponding to soliton solutions) and the number of zeros of the function uj. Both the 
soliton and black hole families of solutions are unstable under linear, spherically symmetric 
perturbations^, with the number of unstable perturbation modes of the solutions being twice 
the number of zeros of of 


Many generalizations of the original spherically symmetric su(2) solitons and black holes 
have been considered in the literature (some of which are reviewed in Ref. Q). For example, 
numerical solutions have been found which retain the spherical symmetry of the original 
solutions but enlarge the gauge group to su(iV) (see, for example. Ref. 0). The solution space 
is more complicated with the larger gauge group, but solutions still exist in discrete families. 
Furthermore, all asymptotically flat, spherically symmetric, soliton and black hole solutions 
with arbitrary gauge group are unstable under linear, spherically symmetric perturbations^. 

The model can also be generalized by considering space-times which are not asymptoti¬ 
cally flat or which have more than four dimensions. In four-dimensional space-time, discrete 
families of spherically symmetric soliton and black hole solutions of su(2) EYM also exist in 
asymptotically de Sitter space-time^i, but, like their asymptotically flat counterparts, they 
are unstable!^. If one considers higher-dimensional space-times, in order to have spherically 
symmetric hnite mass solutions, the YM action must be modihed by the addition of higher- 
order curvature terms^^. With these additional terms, soliton and black hole solutions have 
been found in both asymptotically flat and asymptotically de Sitter space-times^. 

About ten years after the discovery of four-dimensional, spherically symmetric, purely 
magnetic, asymptotically flat, solitons and black holes in 5u(2) EYM, their analogues in 
four-dimensional asymptotically anti-de Sitter (adS) space-time were foundi^^— . The purely 
magnetic su(2) gauge held is still described by a single function u, but now continuous 
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families of solutions are found, which are indexed by the event horizon radius Vh as before 
(including r/j = 0 for solitons), the negative cosmological constant A and the value of the 
gauge field function on the horizon ojh (there is an alternative parameter for soliton solutions, 
which governs the behaviour of the magnetic gauge field function near the origin). One strik¬ 
ing feature of the families of solutions is the existence, for sufficiently large |A|, of solutions 
where the magnetic gauge field function u has no zeros. These solutions where u is nodeless 
are particularly important because at least some of them are stable under linear, spherically 
symmetric, perturbations^^'— . The existence, for sufficiently large |A|, of soliton and black 
hole solutions which are stable under linear, non-spherically symmetric, perturbations has 
also been provert^'^^. In this paper we consider only four-dimensional, spherically symmetric 
solutions, but asymptotically adS generalizations to higher-dimensions^^ or non-spherically 
symmetric space-times^iiSS do exist. 

One generalization which has received a great deal of attention in the literature over 
the past seven years is topological EYM black holes in adS, in particular the relevance of 
black holes with planar event horizons to models of holographic superconductors (see, for 
example, the recent review^ for more details and references). Purely magnetic black holes 
with non-spherical event horizon topology in 5u(2) EYM in adS appeared in the literature 
soon after their spherically symmetric counterparts^^. Unlike the situation in asymptotically 
flat space-time^^, in asymptotically adS space-time su(2) EYM black holes and solitons can 
have nontrivial electric and magnetic fields. While spherically symmetric dyonic solutions 
(both solitons and black holes) in su(2) EYM in adS were found soon after the purely 
magnetic black holesl^, topological dyonic solutions have been studied only more recently. 
Gubser— considered four-dimensional dyonic su(2) EYM black holes in adS with planar 
event horizons. He found a second-order phase transition between the embedded planar 
Reissner-Nordstrom-adS black hole and a black hole with a nontrivial YM field condensate. 
Planar EYM black holes in adS have sub seq uently been widely studied as models of p-wave 

for a selection of work in this area). 


superconductors^! (see also Refs. 


23 and 


28 


Returning to four-dimensional, spherically symmetric, purely magnetic, asymptotically 
adS solutions, a natural question is whether the above stable su(2) solitons and black holes 
have generalizations with a larger 5u{N) gauge group. The answer is affirmative: such 
solutions have been found numerically for gauge groups 5u(3) and su(4)!S. For the larger 
5u{N) gauge group, the purely magnetic gauge field is described by Y — 1 functions Uj 
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(see section Hi Al below). As in the su(2) case, there are continuous families of solutions, 
parameterized by the negative cosmological constant A, the event horizon radius Vh (with 
r/j = 0 for soliton solutions) and A^ — 1 parameters describing the form of the gauge held 


functions either on the event horizon or near the origin. Numerically it is found that, if 


|A| is sufficiently large, then there are solutions in which all the gauge held functions ojj 
have no zeros. For general N, the existence of such nodeless, spherically symmetric, purely 
magnetic, asymptotically adS, soliton and black hole solutions of su(A^) EYM has been 
proven for sufficiently large |Ap^. 

In this paper we address the question of whether these soliton and black hole solutions 
of su(Y) EYM in which all the magnetic gauge held functions have no zeros are stable. The 
outline of the paper is as follows. In section [IT] we introduce su(Y) EYM with a negative 
cosmological constant and the ansatz^ for a spherically symmetric gauge potential. We 
derive the held equations describing static, purely magnetic, conhgurations and the per¬ 
turbation equations for linear, spherically symmetric, perturbations. With an appropriate 
choice of gauge, the perturbation equations decouple into two sectors: the sphaleronic and 
gravitational sectors. These are considered in sections IIIII and HVl respectively. Finally we 
present our conclusions in section |Vl 

II. THE EINSTEIN-YANG-MILLS EQUATIONS 
A. Action, metric and gauge potential 

The action for four-dimensional su{N) Einstein-Yang-Mills (EYM) theory with a negative 
cosmological constant A < 0 is: 



( 2 . 1 ) 


where R is the Ricci scalar, is the non-Abelian gauge held and Tr denotes a Lie algebra 
trace. Throughout this paper, the metric has signature (—, -|-, -|-, -|-) and we use units in 
which AttG = 1 = c. In addition, the gauge coupling constant is hxed to be equal to unity. 
Varying the action fl2.ip yields the held equations: 





( 2 . 2 ) 
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where the Yang-Mills stress-energy tensor is 


Tf,u = Tr 


- -g^^FxaF 


Xa 


(2.3) 


which involves a Lie-algebra trace. The Yang-Mills gange held F^i, is given in terms of the 
gange potential by 

F^i, = dfj^Ay — + [^4^, A^]. (2-4) 

Onr focns in this paper is on eqnilibrinm, static, spherically symmetric, soliton and black 
hole solntions of the held eqnations fl2.2p and time-dependent, spherically symmetric, pertnr- 
bations of those eqnilibrinm solutions. We therefore consider a time-dependent, spherically 
symmetric, geometry, whose metric in standard Schwarzschild-like co-ordinates takes the 
form 


ds^ =dti^ + fi ^ dr"^ + siv? 9 dcjP', (2.5) 


where the metric functions /r(t, r) and S'(t, r) depend on the co-ordinates t and r only. Since 
we have a negative cosmological constant A, it is useful to write the metric function 
in the form 


/i(t,r) = 1 - 


2m{t,r) AF 


( 2 . 6 ) 


r 3 

In our later analysis we will also hnd it useful to dehne another function A{t,r) such that 


S(t, r) = exp A{t, r). 


(2.7) 


With this metric ansatz the relevant components of the Einstein tensor are: 

~ ^ + a) ’ 

Gtr = - —, 

nr 

Grr = in’Sr + 2S'nr - 5 + /i5), (2.8) 

nor^ 

where here and throughout this paper we use a dot to denote d/dt and a prime to denote 
d/dr. Note that we do not need to consider the Gqq or G^^ components of the Einstein 
tensor as the held equations involving these components follow from those involving the 
components in fl2.8p by the Bianchi identities. 

We make the following ansatz for a time-dependent, spherically symmetric, su(Y) gauge 
potential^: 

A = Adt + Bdr + ^{G-G^) d9 - ^ [{G + G^) sm9 + D cos9] (2.9) 
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where A, B, C and D are all {N x N) matrices depending on the co-ordinates {t,r) only 
and is the Hermitian conjngate of C. With this gange potential ansatz, the non-zero 
components of the gange held are: 


Ftr = B- W, 

Fto = l{{C-C^)-+[AC-C^]}, 

= + C^)' + [A, C + C^] } sin e, 

Fre = l{{C-C^y + [B,C-C^]}, 

Pr^ = ~{iC + C^y +[B,C + C^]} sin e, 

Fe<p = ~ { [C*, C^]-D} sin 6. (2.10) 

In compnting the component Fo^p we have made nse of the identities^ 

[D, C] = 2C, [D, C^] = -2C^. (2.11) 

The matrices A and B are diagonal and traceless, and we dehne fnnctions aj{t, r) and r) 
for j = 1,... snch that 


A = i Diag (ai(f, r),... aN^t, r)), 

B = iDiag(/3i(t,r),. ..l3N{t,r )), (2.12) 

where the fact that these two matrices mnst be traceless means that 

N N 

= 0 = (2.13) 

i=i i=i 

The matrix C is npper triangnlar, with non-zero entries only immediately above the main 
diagonal. These entries are given in terms of fnnctions u}j{t, r) and (t, r) for j = 1,, N—1 
by 

C',,,+i = a;,(f,r)e*^^(*’^). (2.14) 

Finally, the matrix H is a constant diagonal matrix—: 

D = Diag (A^ - 1, iV - 3,... , -iV + 3, + 1). (2.15) 
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B. Static solutions 


For static solutions, all field variables depend only on the radial co-ordinate r. We denote 
static equilibrium functions with a bar (e.g. Uj) to distinguish them from the time-dependent 


perturbations which we shall consider shortly. The static equilibrium solutions in which we 
are interested are purely magnetic, which means that we set the electric gauge field functions 
aj{r) = 0 for all j = 1,..., The remaining gauge freedom is then used to set all the 
functions /3j(r) = 0 for j = 1,..., N—. From now on we assume that none of the magnetic 
gauge functions ojj{r) are identically zero. In asymptotically flat space, other families of 
solutions have been found when this assumption is relaxed^. Assuming that none of the 
ojj{r) are identically zero, one of the Yang-Mills equations becomes^: 


7)('') = 0, j = 1,...,N -1, 


(2.16) 


and the gauge held is described by the N — 1 magnetic gauge held functions u!j{r), j = 
1,..., N — 1. We comment that our ansatz fl2.9l) is by no means the only possible choice in 


su{N) EYM (in Ref. [32 all irreducible models are explicitly listed for N < 6, and techniques 


for hnding all spherically symmetric su(Y) gauge potentials are developed). 

1. Static field equations 

For purely magnetic, static equilibrium solutions as described above, the held equations 
fl2.2p simplify as follows. The Einstein equations take the form: 


m' = /i(r)F -|- r^n. A' 


2F 


(2.17) 


r 


where 


N-l 






(2.18) 


The N — 1 Yang-Mills equations take the form 



(2.19) 
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where 


( 2 . 20 ) 


- 1 - + 2 (^i-i + ^i+i) • 

The held equations fl2.17l I2.19P are invariant under the transformation 

u}j{r) -£hj(r) (2.21) 

independently for each j, and also under the substitution: 

j^N-j. (2.22) 


2. Boundary conditions 


The held equations (12.171 [2T^ are singular at the origin r = 0, a black hole event horizon 
r = Th (where fl{rh) = 0) and at inhnity r —)■ cxd. Below we briehy outline the form of the 
equilibrium held functions in a neighbourhood of the singular points. The existence of local 


solutions near these singular points, with the forms below, is proven in Ref. 


30 


a. Origin The form of the static held functions near the origin is rather complicated. 
In particular, to completely specify the gauge held in a neighbourhood of the origin, a 
power series up to 0{r^) is required, involving N — 1 initial parameters. These N — 1 
parameters, together with the cosmological constant A, completely determine the solution 
in a neighbourhood of the ori gin. The details of this power series can be found in Ref. l30 


(following the analysis of Ref. l33| for the asymptotically hat case). For our analysis in this 
paper, we only require the leading order behaviour of the static held functions, which is: 


m(r) = + O(r^), 

S{r) = So + S2r^ + Oir^), 

a;,(r) = ±[j(iV-j)]^ + 0(r2), (2.23) 


where m 3 . So and S 2 are constants. Without loss of generality, we take the positive sign 
in the form of ujj{r) due to the invariance of the held equations under the transformation 

dMH). 

b. Event horizon Assuming there is a non-extremal black hole event horizon at r = rh, 
the metric function Ji{r) will have a single zero there. This hxes the value of m(r/i) to be 

f (2.24) 
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In a neighbourhood of the horizon, the held variables have the form 

fh{r) = fh{rh) + rh'{rh) (r — Vh) + O {r — Thf , 
u}j{r) = u}j{rh) + u'j{rh) (r - tk) + O {r - Vhf , 

S{r) = S{rh) + S'{rh) (r - rh) + O {r - ru ), (2.25) 

where fri'{rh), oj'jirh) and S'{rh) can be written in terms of the constants 0 Jj{rh) and Sir^) 
by using the held equations (12.17112.19p . Again, due to the invariance of the held equations 
under the transformation fl2.2ip . we may take cjj{rh) > 0 without loss of generality. The 
N — 1 initial parameters u}j{rh), together with the cosmological constant A and event horizon 
radius completely determine the solution of the held equations in a neighbourhood of 
the horizon^. 

c. Infinity As r ^ cxo, the held variables have the form: 

fh{r) = M + O , 

S{r) = l + 0{r-^) , 

Uj{r) = Uj^oo + Cjr~^ + O (r~^) , (2.26) 

where M, and Cj are constants. 

3. Embedded solutions 

Despite the complexity of the static held equations (12.17112T^ . there are some embedded 
solutions which will be useful in our later analysis. 
a. Schwarzs child-adS If we set 

u,{r) ^ ±.yj{N-j), (2.27) 

for all j = 1,..., A^ — 1 then the components of the gauge held strength tensor fl2.1UI) vanish 
identically. In this case the stress-energy tensor fl2.3p therefore also vanishes and we obtain 
the Schwarzschild-adS black hole solution with 

m(r) = M, S{r) = 1, (2.28) 

where M is a constant representing the mass of the black hole. Setting M = 0 gives pure 
adS space-time as a solution of the held equations. 
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b. Reissner-Nordstrom-adS Alternatively, if we set 


ujj{r) = 0 


(2.29) 


for all j = 1 ,..., — 1 then the gauge held strength tensor fl2.10l] does not vanish as has 

a contribution from the nonzero matrix D fl2.15p . In this case we obtain the magnetically 
charged Reissner-Nordstrom black hole solution with 


m{r) = M — 


2r 


S{r) ^ 1, 


(2.30) 


where the magnetic charge Q is hxed to be 

= liV (AT + 1) (AT - 1). (2.31) 

0 

c. Embedded su(2) solutions The above two solutions are effectively Abelian embedded 
solutions. For all N > 2, there is another class of embedded solutions, corresponding to 5u(2) 
non-Abelian solutions. To obtain these solutions, we write the N — 1 magnetic gauge held 
functions obj{r) in terms of a single magnetic gauge held function u}{r) as follows: 


uj{r) = ±u{r)^yj {N -j). 


(2.32) 


It is shown in Ref. 


30l that, by a suitable rescaling of the other held variables, in this case 


the static held equations (12.171 I2.19P reduce to those for the su(2) case with uj{r) as the 
single magnetic gauge held function. Therefore any solution of the su(2) held equations can 
be embedded as a solution of the su(A^) held equations. In particular, setting a)(r) = 1 gives 
the Schwarzschild-adS solution of the embedded su(2) held equations and setting 0{r) = 0 
gives the magnetically charged Reissner-Nordstrom-adS black hole. 


4- Non-embedded solutions 

Genuinely su(A^) static soliton and black hole solutions, which do not fall into one of 
the categories described in section III B 31 have been studied in some detail already in the 
literature^^. Therefore in this section we very briehy describe some of the key features of 
the solutions which are required for our subsequent analysis. 

As discussed in the introduction, continuous families of solutions of the held equations 
(12.171 12.191) are found numerically. The solutions are parameterized by the cosmological 
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FIG. 1. Example nodeless solutions for su(3) EYM with A = —10. In each case we plot the gauge 
field functions wi(r) and uj 2 (r) (the typical behaviour of the metric functions can be found in the 
example solutions plotted in Ref. |29|). In (a) we show a soliton solution, and in (b) a black hole 
solution with = 1. 


constant A, the event horizon radius (we can consider r/^ = 0 to represent soliton solutions) 
and, for su(A^), there are a further N — 1 parameters which describe the gauge held (see 
section HlB 2j) . For black holes, these N — 1 parameters are simply the values of the gauge 
held functions on the horizon 0Jj{rh) (I2.25|l . For soliton solutions the si 


29 


nation is more 
and 30 


complicated, details of the parameters in this case can be found in Refs. 

In Ref. [30 the existence of genuinely su(A^) solutions of the static held equations fl2.17l 


I2.19P in a neighbourhood of the above embedded su(2) solutions was proven for all N, for 
sufficiently large |A|. In this article we focus on those su(A^) solutions which are close to 
embedded su(2) solutions and for which all the gauge held functions ujj{r) have no zeros. 


In Refs. 


29 


and 


34l we have presented various phase space plots for 5u(3) and 5u(4) which 


demonstrate numerically the existence of regions of these nodeless soliton and black hole 
solutions. Here we simply plot, in hgures [T] and |2l examples of nodeless soliton and black 


hole solutions for su(3) and su(4) respectively, referring the reader to Ref. 
details of the phase space of solutions. 


29 


for further 


C. Perturbation equations 

In this paper we are interested in linear, spherically symmetric, perturbations of the static 
equilibrium solutions discussed in section IIIBI Our particular interest is in time-periodic, 
bound state, perturbations which vanish at either the origin or event horizon (as applicable. 
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FIG. 2. Example nodeless solutions for 5u(4) EYM with A = —10. In each case we plot the gauge 


field functions a;i(r), W 2 (r) and U 3 {r) (t 
in the example solutions plotted in Ref. 
hole solution with rh = 1. 


le typical behaviour of the metric functions can be found 


29|). In (a) we show a soliton solution, and in (b) a black 


for soliton and black hole solutions respectively) and at infinity. 

To derive the perturbation equations, we write our time-dependent held variables as a 
sum of the static equilibrium quantities (denoted by a bar, e.g. / 2 (r)) plus small perturbations 
(denoted by a 6, e.g. 6fi{t,r)) as follows: 

p(f,r) = / 2 (r) + (5p(f,r), 

S{t,r) = S{r) + 6S{t,r), 
m{t, r) = m(r) -|- 5m(f, r), 

A(f, r) = A(r) -b 5A(f, r), 
aj{t, r) = 6aj(t, r), 

= Sf3j{t,r), 

ujj{t,r) = Oj{r) + 6uj{t,r), 

-fj{t,r) = (2.33) 

Recall from section fll Bl that the gauge held functions aj, jSj and 7 ^ all vanish for static 
equilibrium solutions, but here we consider non-zero perturbations of these parts of the 
gauge potential. 

The perturbation equations are found by substituting the held variables in the form fl2.33p 
into the components of the Einstein tensor fl2.8p and gauge held fl2.10l) . and then working 
out the held equations fl2.2p . We work only to hrst order in the perturbations and simplify 
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the resulting equations using the static equilibrium held equations fl2.17l I2.19p . 


First of all, the linearized Einstein perturbation equations become 


5/i' 

Sfi 

6A' 


[-5/i - 2T6n - 2j26T - 2r‘^6U] , 
r 

Ur, 


(2.34a) 

(2.34b) 

(2.34c) 


where 


N-l 

sr = 2j2'^'M. 

i=i 
1 ^ 

5n = — ^ [ujj — - 1 + 2j] [ujjSuj - Uj-x5uj-i], 

i=i 

Af-l 

SH (2.35) 

j=i 

and we remind the reader that we are using a dot to denote d/dt and a prime to denote 
d/dr. 

Given the form of the elements of the matrix C fl2.14l) . it is useful to consider the following 
combinations of the perturbations 6uj(t, r) and 57j(t, r), for j = 1 ,... — 1: 


S'dj{t,r) = 6uj{t,r) + iUj{r)S'yj{t,r), 

5ip*{t,r) = Suj{t,r) - iu)j{r)S'jj{t,r), (2.36) 


in terms of which the entries of the matrix C are, to hrst order in the perturbations: 


Cj,j+i = ujj{r) + 6^/j{t,r). 


(2.37) 


In terms of 5'ipj-, Si/j*, the linearized Yang-Mills perturbation equations are: 
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0 = /i 


5/3' — 6a" + (sPj — 6a j 


2 _ S' 
r~~S 


+ 


2 r 2 L 


Uj ( 6ipj — 6ip' 


-uij-i ( 6'ijjj_i - 6i)j_i ) + 2cj {6aj - 6aj+i) - 2u -^ {6aj_i - 6aj) 


(2.38a) 


0 = 


fiS^ 




Ui 




w,--i 


U - ^r--) - a', (H’, - H’;) 


0 = 


{6'ipj-i — 6'i6*_i) + 2ci;|_i {6Sj — 6Sj-i) — 2,u‘j (5/3j+i — 5/3^)] , 
6'ipj + ojj {6aj — 5dj+i) + —a)''5/i + ^ ~ 


(2.38b) 


2 / 2^2 L 
+2a}' (5/3j — 5/3j+i)] + 


wl- 


5/i' + /i5 


2 




+ ^ [#' + Qj (ift - 


0 = 


1 

+W(,-5V',], 

1 


{6i)j + 5V’i) + (^^i+i + 5^/’*+i) + 


(2.38c) 


2 /i ^2 L 

+2a;' (5/3j — 5/3j+i)] 


-5^; + cu,- {6a, - 5a,+i) - ic5''5/i + | [-5^/' + cb, (5/3' - 5/3'+,) 




5/i' + /i5 






s / + 7'^^ 


+ 


2 L 


-63p*- + w,- (5/3,- - 5/3,+ 1 ) 


1 

-w^,5^;]. 


a;| (5^5i + 5^5*) - {6'tpj+i + 5V’;+i) - (5?5j-i + 5^5,*_i) 


(2.38d) 


The equations (I2.38al - [2.38bll come from the t and r Yang-Mills equations, respectively, and 
there are N of each of these equations, corresponding to j = 1,... N. The equations (I2.38ch 
I2.38dl) come from taking the real and imaginary parts of the 6 Yang-Mills equation (the 0 
Yang-Mills equation gives the same pair of equations), assuming that all perturbations are 
real, and there are iV — 1 of each of these equations, corresponding to j = 1,... Y — 1. 

Our time-dependent, spherically symmetric gauge held ansatz fl2.9p has a residual gauge 
degree of freedom. For a diagonal matrix g(t, r), consider the following gauge transformation: 


A ^ A + Q ^ 0 , 

B ^ B + Q ^ 0 ', 

C-C^ ^ 0 -^ [C - C^) 0 , 

C + C^ ^ 0 -^ {C + C^) 0 , (2.39) 

under which the gauge held transforms as 

-2 Q~^Ff,^Q. (2.40) 
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We choose the diagonal matrix g so that ^ + g = 0, which enables us to set the pertur¬ 
bations 6aj{t, r) = 0 for all j = 1,. .. N. 

With this choice of gauge, the perturbation equations fl2.34af[^34cll2.38al - l2.38dj) decouple 
into two sectors. The hrst sector contains the Yang-Mills perturbations S^j, j = 1,... ,N and 
j = 1,..., iV — 1 and does not contain any metric perturbations. This sector is known as 
the sphaleronic sector-. This terminology arises from the fact that the su(2) EYM solitons^ 
and black holes^ in asymptotically flat space possess instabilities in this sector->^>^ analogous 
to the unstable mode of the Yang-Mills-Higgs sphaleron^^. The second sector contains 
the perturbations of the metric functions 6fi and 6A and the Yang-Mills perturbations 
6uj, j = 1,..., N — 1. This sector is known as the gravitational sector. As the static 
equilibrium solutions are purely magnetic and spherically symmetric, they are invariant 
under a parity transformation. As a result of this additional symmetry, the two decoupled 
sectors of perturbations transform in a particular way under a parity transformation: the 
perturbations in the sphaleronic sector have odd parity and change sign under a parity 
transformation; the perturbations in the gravitational sector have even parity and do not 
change under a parity transformation. 


In the analysis of the sphaleronic and gravitational perturbation sectors in sections IIIII 
and [IV] respectively, we will change our independent radial variable to the usual ‘tortoise’ 
co-ordinate r*, dehned by 


dr* _ 1 
dr p,S 


(2.41) 


For perturbations of static soliton solutions, we choose the constant of integration such that 
r* = 0 at the origin where r = 0. In this case r* has a maximum value, Tc, as r —)■ cxo. For 
perturbations of static black hole solutions, we choose the constant of integration such that 
r* —)■ 0 as r —)■ oo, and then r* —)• —cxo as the event horizon is approached, r rh- 


III. SPHALERONIC SECTOR PERTURBATIONS 


The sphaleronic sector consists of the odd parity Yang-Mills perturbations djdj (j = 
1,... ,N) and (j = 1,..., A — 1). In the gauge Saj = 0, j = 1,..., N, the sphaleronic 
sector perturbation equations are fl2.38al I2.38bj) and a third perturbation equation which 
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comes from adding equations fl2.38cp and fl2.38dp . The equations are: 


0 = /i 






(3.1a) 


0 = ^ - 7-1 Ws - + 7 ('57+1 - '5ft)]. 




‘ZfJjOOj -\~ fi UJj IXLiJj 


J 


[d'y'j + 6(3j — S(3j+i\ 


(3.1b) 


(3.1c) 


A. Sphaleronic sector perturbation equations in matrix form 


We first simplify these equations by changing our radial co-ordinate from r to the tortoise 
co-ordinate r* fl2.4ip and by introducing new dependent variables Sej (j = 1,...,A) and 
5«h,- (j = l,...,A-l)by: 


5ej = ry/Jl6Pj, 6^j = oJjS'jj. 
The perturbation equations fl3.1af[3TcP then take the form 
0 = 5kj -\- 


6is>, = aj 


(3.2) 


■ ^ (ft'57 - ■^i-fth-i). 

(3.3a) 

+ idr,Qj) 5^j - {dr,UJj-l) 5<hj-l] 


5ej) - uj]_^ - (5ej_i)] , 

(3.3b) 


r* ^ 


UJj 


“1“ hujjdf^ i.^^j — T {h^j) T ((fcj — > 


(3.3c) 


where we have introduced the quantity 

h = (3.4) 

r 

We now express the perturbation equations in matrix form by dehning 

Se = {5ei,...,6eN)^, 

= (3.5) 
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in terms of which the perturbation equations fl3.3af[3y^ take the form 


(3.6a) 

(3.6b) 

(3.6c) 


0 = dr, {h-^5k) + 

6i = h^lC6e - h [7dr,6^ - {dr,J^) , 

= dlS^ + hj^^drje + XSe + W6^, 


where we have dehned an iV x [N — 1) matrix an x iV matrix /C, an [N — 1) x [N — 1) 
matrix W and an {N — 1) x N matrix X as follows: 




^ bJ\ 0 0 ■ ■ ■ 

—ijJ\ U2 0 ■ ■ ■ 

0 — UJ 2 OJ 3 ■ ■ ■ 


0 ^ 

0 

0 


U)N -1 


K, = 


y 0 000 —cujv—1 / 

^ —ujI uj\ 0 • • ■ 0 ^ 

bj\ —bj\ — bj\ ■ ■ ■ 0 

0 bj\ —ijj\ — oj\ • • ■ 0 

: : : ■. ,-,2 

y 0 0 0 ^N-i ~^N-i / 


W = h‘^ Dmg{Wi,...,WN-i), 

A” = 2hdr,J=''^ + {dr,h) 


(3.7a) 


(3.7b) 

(3.7c) 

(3.7d) 


where the quantities Wj, j = 1,..., — 1 are given by (12.201) . Finally, we introduce a vector 

of dimension 277 — 1 by 

in terms of which the hrst perturbation equation (13.bah takes the form 


= dr. 


1 

( In 

o\ . 



\ » 

V J 



and the remaining equations fl3.6bf[3l6c|) can be compactly written as 


(3.9) 


- i =W4', 


(3.10) 
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where we have dehned the operator 


= 



hdr^J= 

X W 


(3.11) 


and X„ denotes the n x n identity matrix. 

It is straightforward to show that the operator U fl3.11l) is real and symmetric when acting 
on pertnrbations which vanish at either the origin or event horizon (as applicable) and at 


inhnity. However, as noted in Ref. the operator U is not elliptic and so the pertnrbation 
eqnation (Id.lOj) is not currently in hyperbolic form. 

For time-periodic perturbations for which ^{t,r) = e*°'*^'(r), the perturbation equations 
(I3.in|l take the form 

(3.12) 


If we can show that the operator W is a positive operator, then the eigenvalues must also 
be positive. This means that a is real and the perturbations are periodic in time. In this 
case small perturbations remain small and there are no unstable modes in the sphaleronic 
sector. Our aim for the remainder of this section will be to show that there are at least some 
equilibrium su(A^) soliton and black hole solutions for which W is a positive operator. 


B. The Gauss constraint 


The hrst of the linearized Yang-Mills perturbation equations fl3.9l) is known as the Gauss 
constraint. A lengthy calculation reveals that the Gauss constraint propagates, in other 
words the perturbation equations (l3.6bf[T6c|l imply that 

= 0 (3.13) 

independently of the Gauss constraint. Equivalently, we may writel^ 


gu = o. 


(3.14) 


Following Ref. Il0|, we integrate (13.hh with respect to time, and choose the constant of inte¬ 
gration (in this case an arbitrary function of r) so that: 




(3.15) 
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which we will call the strong Gauss constraint. Suppose we have a vector of perturbations 
which satisfy the strong Gauss constraint at initial time t = 0, and which initially satisfy 
the Gauss constraint fl3.9p . By virtue of fl3.14p . this vector of perturbations will satisfy the 
strong Gauss constraint at all subsequent times. 

Now consider any vector of perturbations ^ (satisfying the perturbation equations) and 
write it as the sum of two parts: the first, satisfying the strong Gauss constraint and 
the second, ^ 2 , failing to satisfy the strong Gauss constraint. It is shown in Ref. ^ that 
the second vector of perturbations, ^^ 2 , is pure gauge, having the form 


^^2 = 


(3.16) 


where 


, , , 0 \ / 0 0 

0 0 / \ 0 


(3.17) 


is the adjoint of the operator Q fl3.9p . Such perturbations correspond to infinitesimal gauge 
transformations of the form fl2.39p with (for small e) 


0 = exp ( —eT 


(3.18) 


where T is an iV x iV matrix of the form 


t = Diag(Ti,...,T„) 


(3.19) 


with Tj, j = l,...,iV the first N elements in T. Therefore a vector of perturbations 
T' which satisfy the Gauss constraint but not the strong Gauss constraint can be gauge- 
transformed to a vector of perturbations satisfying the strong Gauss constraint. Without 
loss of generality, we may therefore restrict attention to physical perturbations satisfying 
the strong Gauss constraint, which is essentially an initial condition. Since this is a gauge 
transformation of initial data only, the matrix g fl3.18p is time-independent, and so, by 
(I2.39p . this transformation preserves the gauge condition 6aj = 0. 


C. An alternative form of the operator U 

In order to prove the existence of static solutions which have no unstable modes in the 
sphaleronic sector governed by the equations fl3.10p . in the next subsection we will want to 
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show that the symmetric operator lA fl3.11l) is positive. In this subsection we will write the 
operator in an alternative form which will enable us to hnd static equilibrium solutions for 
which W is a positive operator. 


In particular, following Ref. 


10 


, we seek operators y and V such that we may write 


u = x\ + v-g^h^g, 


(3.20) 


where the operator g is given by fl3.9p . We then hnd 


g^h^g 


( Xn Q 
0 0 




h-^dlh -drAhJ^)-{drMJ^ \ 
-{dr^h)J^ y 

(3.21) 


Next dehne the operator y as 

X = dr.X hz, (3.22) 

where Z is some {2N — 1) x (2iV — 1) matrix which does not contain any derivative operators 
and which is to be determined. Then 


= -dl +h{z^- Z) dr, + [h^Z^Z - dr, {hZ)] . (3.23) 


From now on we assume that Z is symmetric, so that there is no hrst order derivative 
operator in yly. Writing the matrix Z in the form 


( Xii Zi2 

y X'12 X'22 


(3.24) 


where Zn is a symmetric NxN matrix, Z 12 is an x {N — 1) matrix and Z 22 is a symmetric 
(A^ — 1) X (A^ — 1) matrix, using (13.21113.23P we hnd that the matrix V dehned in fl3.20p has 
the form 


V 


Vn V 12 

V 21 V 22 




(3.25) 


where 


Vn = -h^JC + {Zf,Zu + ^12^1^2) " dr, {hZn) + h"' (dlh) Xat, 

V 12 = -2hdr,J^ + (XfiZi2 + X 12 XJ 2 ) - dr, (hZn) - 2 {dr,h) 7, 

V 2 I = -2 {dr,K) 7^ - 2hdr,7^ + {Zl,Zu + Z^,Z^^) - dr, {hZ^2) , 

V 22 = -W + {Z^2Xi2 + ^ 2 ^ 2 ^ 22 ) - dr, {hZ22) + 7^7. (3.26) 
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We are free to choose the matrix Z so as to simplify the form of V. We hrst make the choices 


Z\\ — h {dr,h)XN, Z 22 — 0. 
In this case the form of V 12 simplihes to 

V 12 = -2dr, {hj^) - hdr^Zu, 


(3.27) 


(3.28) 


which vanishes if we choose Z 12 such that 

j-rt 

Zi 2 = -2 / h-^dr^ {hJ^) dn, (3.29) 

where r* min = 0 for equilibrium static soliton solutions and r* min = —00 for equilibrium 
static black hole solutions. With this choice of Z 12 , it is straightforward to see that V 21 also 
vanishes. 

The matrix V is then block diagonal, with its diagonal entries being 

Vn = -h^lC + h^Z^2Zl^ + {h-^dr,hYXN, 

V 22 = -W + h^Zl^Zi2 + (3.30) 


D. Conditions for no instabilities in the sphaleronic sector 

For physical perturbations satisfying the strong Gauss constraint fl3.15p . the operator U 
fl3.20p appearing in the sphaleronic sector perturbation equations fl3.12p reduces to 

U = xh + V. (3.31) 

Since the matrix V is symmetric, the operator lA is symmetric and real. From the form of 
the operator x'x 113.23p . the operator fl3.3ip is elliptic. Furthermore, since x^X is a positive 
operator, to show that W is a positive operator it suffices to show that V is a positive matrix. 
The matrix V is block diagonal, and hence positive if its two non-zero diagonal blocks Vn 
and V 22 113.30p are positive. 

Let us begin with Vn. The second and third terms in Vn are manifestly positive, so it 
remains to consider the term —where the matrix K is given by fl3.7bp . For an arbitrary 
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vector X = (xi,..., x^)'^, we have 

- x'^JCx = UjIxI + (ul + U 2 ) xl + . . . + {u %_2 + (^N-i) x\-i + ^%-ix]^ 

— 2(jj‘lxiX2 — 2 uJ2X2X^ — ... — 2a;^_j^XAr_iXAr 
= u\ (xi - X 2 f + ^2 (^2 - Xsf + . . . + (xat.i - Xjvf 

> 0. (3.32) 

Therefore Vn is positive. 

For V 22 , again the second and third terms are manifestly positive. The hrst term, —W 
fl3.7cll is a diagonal matrix, which will be positive if and only if its entries are positive. For 
this to be the case, we reqnire Wj < 0 for j = 1,..., iV — 1, where the qnantities Wj are 
defined in fl2.20p . This gives the following set of ineqnalities to be satisfied by the static 
eqnilibrinm solntions for all r: 

oj\>l + -ul, 

0)2 > 1 + - + qI) , 

- ^ + 2 ’ 

^"n-i > 1 + 2^n-2- (3.33) 

If the ineqnalities fl3.33p are satisfied for all r, then we can dednce that V 22 is positive. There¬ 
fore the matrix V is positive, and hence the operator W is a positive operator. We deduce 
that physical solutions of the sphaleronic sector perturbation equations (I3.12p must have 
positive, so a is real and the perturbations are periodic in time. Therefore small pertur¬ 
bations remain small and the equilibrium solutions have no instabilities in the sphaleronic 
sector. 

We emphasize that the inequalities fl3.33p are sufficient for an equilibrium solution to have 
no unstable modes in the sphaleronic sector; we have no expectation that these inequalities 
are necessary for stability. Our interest in this paper is in proving the existence of stable 
soliton and black hole solutions of the 5u{N) EYM equations. Therefore we will have 
achieved this aim, at least for the sphaleronic sector, if we can find equilibrium solutions 
satisfying fl3.33p for all r. 
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E. Special cases 


Before proving the existence of non-trivial su(A^) equilibrium solutions which have no 
instabilities in the sphaleronic sector, in this subsection we consider the embedded solutions 
discussed in section III B 31 


1. Schwarzschild-adS 

Setting Uj = j (iV — j) and m{r) = M, we hnd that Wj = 0 for all j = 1,..., — 1 and 

so the matrix W vanishes identically. In this case the matrix V is manifestly positive and the 
operator lA is positive when acting on physical perturbations. Therefore the Schwarzschild- 
adS solution, as expected, has no instabilities in the sphaleronic sector. 


2. Reissner-Nordstrom-adS 

If Uj = 0 for all j = 1,..., — 1, then the sphaleronic sector perturbation equations 

reduce greatly. The perturbations 5<I)j 03.21) vanish identically, leaving only the Sej pertur¬ 
bations. The only 6ej perturbations which then satisfy the strong Gauss constraint have the 
form 

6ej = Zjh = Zj ^ (3.34) 

where Zj are arbitrary constants. These perturbations do not vanish at the origin or the 
event horizon or inhnity unless Zj = 0 for all j = 1,..., A^. For bound state perturbations 
which vanish at either the origin or the event horizon (as applicable) and at infinity, the 
only possibility is Zj = 0 and hence 5ej = 0. This means that there is no dynamics in the 
sphaleronic sector when the static equilibrium solution is embedded, magnetically-charged. 
Abelian Reissner-Nordstrom-adS. The only allowed perturbations of the gauge potential 
correspond to gauge transformations. 
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3. Embedded 5 u( 2 ) solutions 


Setting ujj = uJ{r)^yj {N — j), the sphaleronic sector perturbation equations simplify 
considerably. In particular, the matrix W fl3.7cp reduces to 

W = h^l-od^)lN_i. (3.35) 

The inequalities (I3.33p are then all satished if 

u{rf > 1 (3.36) 


for all r, in which case the operator It is positive and there are no instabilities in the 
sphaleronic sector. 

It remains to prove the existence of embedded su(2) solutions satisfying fl3.36p for all r. 
Embedded su(2) black holes are parameterized by the radius of the event horizon the 
cosmological constant A and a)(r/j). Fix and choose oj{rh) > 1. From the Yang-Mills 
equation (12.1911 in the su(2) case, we have oj'irh) > 0, so that oj is an increasing function of 
r close to the horizon. Embedded su(2) solitons are described by a single parameter 6 , such 
that, near the origin, 

a;(r) = 1 -I--I- O(r^), (3.37) 


where, without loss of generality, we are assuming that a}(0) > 0. Choose 6 > 0, so that 
a;(r) > 1 in a neighbourhood of the origin and > 0 for r sufficiently small. Also from 
fl2.19p . we see that the gauge function u cannot have a maximum if a) > 1. Therefore u 
will be an increasing function of r for all r > Vh for our black hole solution and all r > 0 
for the soliton solution. Therefore fl3.36p will be satished and these solutions will have no 
instabilities in the sphaleronic sector. The same argument applies if a) < — 1 either at the 
horizon or near the origin. 

n 

In Ref. Il5|it is proven that 5 u( 2 ) black holes have no instabilities in the sphaleronic sector 
of su(2) EYM perturbations as long as the g aug e function a;(r) has no zeros, with no further 


conditions on a;(r). While the proof in Ref. 


15 


is for black holes only, the argument carries 


over trivially to the soliton case. We note that here, we have a stronger sufficient condition 
fl3.36p for 5 u(2 ) solutions embedded in su(Y) EYM to have no instabilities in the su(Y) 
sphaleronic sector. This is to be expected since the su{N) EYM sphaleronic sector has more 
degrees of freedom {2N — 2, comprising N functions (3j whose sum must vanish and N — 1 
functions jj) than the su(2) EYM sphaleronic sector (which has just two). 
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F. Existence of static solutions with no sphaleronic sector instabilities 


We now turn to proving the existence of non-trivial su(iV) EYM solitons and black 
holes having no instabilities in the sphaleronic sector. From the above analysis, all that is 
required is to show the existence of equilibrium solutions satisfying the inequalities fl3.33p . 
The argument for both soliton and black hole solutions is straightforward, based on results 
from Ref. 30. 


Due to the symmetry of the held equations fl2.17l I2.19P under the transformation fl2.2ip 
it is sufficient to consider gauge held functions such that uj > 0 near either the origin if 
we are considering a soliton solution or the event horizon if we are considering a black hole 
solution. First we dehne the open region TZ which is the set of all positive values of the 
equilibrium gauge held functions Uj > 0, j = 1,..., N — 1 such that Wj < 0 (so that 
the inequalities fl3.33p are strictly satished for all points in TZ). From the argument in the 
previous subsection, there are embedded su(2) soliton and black hole solutions such t 


tia' 


the values of the gauge held functions he in TZ for all r. From Proposition 9 of Ref. 


30 


there are genuinely (that is, non-embedded) su{N) soliton and black hole solutions whose 
initial parameters (either near the origin or the event horizon, see section III B 21 for details) 


he in a neighbourhood of the initia 
holes. Propositions 3 and 6 of Ref. 


parameters for the embedded su(2) solitons and black 


30 


tell us that the equilibrium gauge held functions ojj 


are analytic functions of the initial parameters and the radial co-ordinate. Fix Vh for the 
black hole solutions under consideration and set ri S> max{l,r/i} (with Vh = 0 for soliton 
solutions). Then, by analyticity, providing our su{N) solutions have initial parameters 
sufficiently close to the initial parameters for the embedded su(2) solutions, the gauge held 
functions ojj{r) for the 5u{N) solutions will remain close to the embedded su(2) solutions 
for all r < ri and hence also within the region TZ for all r < ri. Providing we have chosen ri 
sufficiently large, for r > ri we are in the asymptotic large r regime discussed in section 4.2 


of Ref. 


30- The upshot of that analysis is that, by taking ri sufficiently large, the change in 


the gauge held functions as r —)■ cxo from r = ri can be made arbitrarily small. Therefore, 
since TZ is an open region, the gauge held functions u)j for our su(iV) solutions will remain 
inside TZ for all r > ri. 

By way of illustration, in hgures [3] and IH we show how —Wj fl2.20p depend on r for the 
example soliton and black hole solutions plotted in hgures [1] and [2] respectively. In all cases 
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FIG. 3. —Wi{r) and —W 2 {r) for the example nodeless 5u(3) solutions shown in figuredl (a) soliton 
and (b) black hole. In both cases —Wj > 0 for all r, so that the inequalities ()3.33p are satisfied 
and the solutions have no instabilities in the sphaleronic sector. 



FIG. 4. —lyi(r), —W 2 {r) and —W^{r) for the example nodeless su(4) solutions shown in hgure[2l 
(a) soliton and (b) black hole. In both cases —Wj > 0 for all r, so that the inequalities (|3.33l) are 
satisfied and the solutions have no instabilities in the sphaleronic sector. 

we see that —Wj > 0 for all r, so that these example solutions have no instabilities in the 
sphaleronic sector. 


IV. GRAVITATIONAL SECTOR PERTURBATIONS 

The gravitational sector consists of the even-parity Yang-Mills perturbations duj {j = 
1,... ,N — 1) and the metric perturbations 6fi and 6 A. The governing equations are the 
perturbed Einstein equations fl2.34afl2.34cp together with the linearized Yang-Mills equations 
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formed by subtracting equations fl2.38cp and fl2.38dD . namely: 


5' 


0 = —6ojj + + /rS'^ [^1^' + ) Suj + 5^' + ^5 

flS^ 


S'\ , S' 


|^Pbj(5cUj — 2,ujj6ujj ujj^iU!j6u!j^i “b ujj—iujj5ujj—i^ 


(4.1) 


A. Metric perturbations 

The linearized Einstein equations (I2.34af[2.34cp can be used to eliminate the metric per¬ 
turbations from the remaining gravitational sector perturbation equations fl4.1l) . The lin¬ 
earized Einstein equation fl2.34bp can be immediately integrated to give 

A- 

6fi = —^ ^ 0 )'-b 5E(r), (4.2) 

i=i 

where 6Y (r) is an arbitrary function of r alone. Using the linearized Einstein equation 
(I2.34al) gives, after some lengthy algebra 

6Y' = -1 (1 + 2f) SY, (4.3) 

where T is given by (I2.18p . Integrating (14.31) we find 

6Y{r) = To exp J ^ [l + 2f (r')] (4.4) 

where Tq is a constant. The lower limit on the integral rg = 0 if we are considering pertur¬ 
bations of a static soliton solution, vq = if we are considering perturbations of a static 
black hole solution. Since we require our perturbations to vanish at either the origin or black 
hole event horizon, as relevant, it must be the case that To = 0 and hence 6Y{r) = 0. We 
therefore have 

A- 

611 = —(4.5) 

j=i 

B. Gravitational sector perturbation equations in matrix form 

Now that we have the form (14. 5 h of the metric perturbation dp, together with (I2.34cp 
for the perturbation dA', we can eliminate the metric perturbations from the gravitational 
sector perturbation equations fl4.ip . 
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First we consider the quantity 



where hll is given in fl2.35p . The right-hand-side of fl4.6p depends only on the perturbations 


6ujj and not on their derivatives. Next we dehne a vector of perturbations as follows: 



( 4 . 7 ) 


Changing the radial co-ordinate to the tortoise co-ordinate (I2.4ip . the gravitational sector 
perturbation equations fl4.ip take the form 


— dCj = —dl duj + MSuj. 

r * 


( 4 . 8 ) 


The (iV — 1) X (iV — 1) matrix Ai depends only on the static equilibrium solutions and does 
not contain any derivative operators. To simplify the entries of A4, we make extensive use of 
the static equilibrium held equations (12.17112.19p . After a lengthy calculation, we can write 
the entries of the symmetric matrix A4 as follows. There are three different types of entry 
which have different forms: (i) the diagonal entries Aijj, (ii) entries immediately above and 
below the diagonal Aijj+i and (hi) other entries not on the diagonal nor immediately above 
or below it Aij,k {k ^ j,j + 1). We give these entries explicitly below, where there is no 
summation: 



( 4 . 9 ) 


where we have dehned 


Q = A.p. + isr.s + 

fi S r 


( 4 . 10 ) 

and then 


We now consider time-periodic perturbations for which Suj{t,r) 
the gravitational sector perturbation equations (14.81) are: 


a^Su: = —dl 6u) + 
/ * 


( 4 . 11 ) 


Since (14.111) has the form of a standard Schrodinger-like equation, the operator on the right- 
hand-side of (14.lip is positive if the matrix A4. is positive. If this is the case, then cr^ is real 
and there are no unstable modes in the gravitational sector. 






















C. Special cases 


The gravitational sector perturbation equations fl4.8p are rather complicated in general, 
so first we consider some special cases. 


Schwarzs child- adS 


Setting ojj = {N — j) for j = 1,... — 1, the entries in the matrix Ai in the 

gravitational sector perturbation equations reduce to 








-j(iv-j), 


■^j,k — 0, k ^ j, j 


(4.12) 


Therefore we have 


Ai — 


where Sn-i is the constant matrix with entries 


^N—l,j,k (-^ 'j} ^ ^j+^tk —l,fc] 


It is shown in Ref. 


29 


(4.13) 


(4.14) 


that the eigenvalues of the matrix are k{k + 1) for k = 1,... 1, 


so that the matrix (and therefore the matrix Ai) is positive. The net result of this 

is that, as anticipated, the embedded Schwarzschild-adS solution has no instabilities in the 
gravitational sector. 


2. Reis sner-Nordstrom-adS 

In this case uj = 0 for all j = 1,..., iV — 1, and the matrix Ai reduces to 

M = (4.15) 

Therefore the matrix Ai is negative dehnite everywhere and, as expected, the embedded 
magnetically charged Reissner-Nordstrom-adS solution is unstable. 
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3. Embedded 5 u( 2 ) solutions 


We now have ujj = for j = — 1, and the entries of the matrix 

A4 fl4.9p take the form 




j2S^ 


2fiS^ 


{uj^ -1) +j{N - j) 

I 

= -VjiN-j)ij + l) (N-j-l) 

4Q 


= -^j{N-j)k{N-k) 


fiS'^ _2 , 4Q _ 2 , 8S ^ _ 2 ^ -a 

8 ^ 


LJ 


{dr.Ojf + — (l - 


/iS'r ' ' 

In this case it is helpful to consider the matrix W1 as a sum of three parts: 

A4 = A/i + A /2 + A /3 5 


(4.16) 


(4.17) 


where 




A/3 — 


4Q _.2 S'S' / _2\ -o - 

(dr^o;)-^ (1 - cu j uOr^uj 


fiSr 


^N-l, 


(4.18) 


where the constant matrix is given in fl4.14p and the constant matrix has entries 


SN-i,j,k = ^/j{N-j)k{N-k). 


(4.19) 


The hrst matrix, A/), is positive if a)(r)^ > 1 for all r, which is the same sufficient condition 
as we found previously for no instabilities in the sphaleronic sector (see section IIIIESp . 
Since the matrix Tat-i is positive, the second matrix A /2 is also positive. The positivity of 
the third matrix, A /3 is less clear-cut. However, it has been shown^^ that \/\K\(jj' —)■ 0 for 
all r as |A| —)■ cx3. Therefore, for sufficiently large |A|, the third matrix A /3 is negligible 
compared with A/) and M 2 - Hence, for sufficiently large |A|, if £h(r)^ > 1 for all r, the matrix 
A4 is positive and there are no instabilities for embedded su( 2 ) solutions in the gravitational 
sector of 5\x{N) EYM perturbations. 

it shown that a)(r)^ > 1/3 is a sufficient condition for 5 u( 2 ) black holes to 


In Ref. 


15 


have no instabilities in the gravitational sector of 5 u( 2 ) EYM perturbations provided |A| is 
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sufficiently large. As with the sphaleronic sector perturbations, for embedded su(2) solutions 
in su(A^) EYM, we have a stronger sufficient condition for the absence of unstable modes. 
This is to be expected because of the greater number of degrees of freedom in the su(Y) 
gravitational sector perturbations than in the 5u(2) gravitational sector perturbations. We 
emphasize that the condition of (h(r)^ > 1 for all r and sufficiently large |A| is a sufficient 
condition, and there may be solutions which do not satisfy this condition but which are 
nonetheless stable. 


D. An alternative form of the gravitational sector perturbation equations 


In order to show, in the next section, that there exist both soliton and black hole non- 
embedded solutions of the static equilibrium held equations which have no instabilities in 
the gravitational sector, we shall follow the method of Ref. Il9| and employ a nodal theorem 
due to Amann and Quittner— which allows one to count the number of bound states of a 
Schrodinger-like equation. In this subsection we state Ama nn and Quittner’s result and cast 
our gravitational sector perturbation equations fid.lip in the form required for the application 
of the theorem in section IIVEI We will need to consider solitons and black holes separately. 


1. The nodal theorem 

Amann and Quittner’s theoren>2I is concerned with the number of bound states of a 
radial Schrodinger-like operator. Let D be the linear differential operator 


'Du = — 


d 

dp 


dp 


P" 


-Mp) + ^{P) 


u. 


(4.20) 


acting on n-dimensional vectors u{p), where p G [0, cxo) lies on the half-line. The n x n 
matrices 2l(p), f8(p) and (L(p) are assumed to be real, symmetric, smooth and uniformly 
bounded on [0, oo). It is further assumed that 2l(p) is uniformly positive dehnite on [0, cxd), 
that is, there is a constant o > 0 such that 


2 I(p) > a > 0 for 0 < p < oo, 

and that 23(0) is non-negative. 

The theorem is concerned with the eigenvalue problem 


(4.21) 


Du = Xu, tt(0) = 0, 


(4.22) 
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where u & L 2 ((0, 00 ), M"'). Following Ref. l37|, we further assume that the bottom of the 
essential spectrum of 2) is positive and that the eigenvalue problem fl4.22p has only hnitely 
many negative eigenvalues. Sufficient conditions for this assumption to be valid are^: 


C(p) —0 as p —)• CX) 


(4.23) 


and 

£(p)p2 + ®(p) > -b (4.24) 

for some b < 1/4 and all sufficiently large p. 

The statement of Amann and Quittner’s theorem involves an auxiliary problem, which 
we now state. Choose n linearly independent real, constant, n-dimensional vectors e^-, 
j = 1,..., n. Let c > 0 and let 

if, = [ni,...,nQ (4.25) 


be the n x n matrix whose columns are the solutions of the n initial value problems 


d 

Duj = 0, c < p < 00 , Uj[c) = 0, —Uj{c) = Cj, 

(JjU 


(4.26) 


for j = 1,... ,n. We then define a scalar function ^{p) by 


dip) = detilc(p). 


(4.27) 


We are now in a position to quote Amann and Quittner’s theorem^l; 


Theorem 1 If c> 0 is sufficiently small and d > c is sufficiently large, the number of zeros 
(counted with multiplicities) in the interval (c,b) of the function dip) equals the number of 
negative eigenvalues of U.2S\ ) (counted with multiplicities). 


In order to apply this theorem, we need to cast the gravitational sector perturbation 
equations fl4.1ip in the form fl4.22p . choose a suitable co-ordinate p and check that the 
matrices 2l(p), 23 (p) and (i(p) satisfy the required conditions (14.211 I4.24p together with 


the requirement that iB(0) is non-negative. To do this, as in Ref. 


19 


, we need to consider 


soliton and black hole solutions separately. We consider black holes first as this case is 
simpler. 
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2. Black holes 


Following Ref. 


19 


for static black hole solutions we take 


p = -r* e [0, oo) 


(4.28) 


where r* is the tortoise co-ordinate fl2.4ip . This means that p —)■ 0 corresponds to r —)■ oo, 
and p —)■ oo corresponds to r ^ approaching the event horizon. 

We choose the (iV — 1) x {N — 1) matrices appearing in the differential operator T) fl4.20p 
as follows: 

2t(p)=Xjv-i, ®(p) = 0, (t{p)=M, (4.29) 

where the matrix A4 has entries fl4.9p . Taking a = | > 0, the condition fl4.2ip is automat¬ 
ically satished, and, furthermore, iB(0) =0 is non-negative. It remains therefore to check 
the conditions on the matrix Cf fl4.23f[T?24p . 

We hrst examine the behaviour of the matrix £(p) as p —?• 0, that is, r —)■ oo. Using the 
boundary conditions fl2.26p . and noting that, as r —)■ oo, 

A.C ■ 

dr.i^j = -^ + 0{r~^), Q = -Ar + 0{1), (4.30) 


we hnd that the leading-order behaviour of € is given by the entries 




^ + 9 (^j+l,oo A- ^j-l^oo) 


-I- 0{r 


-I'l 


A 


2 i^yoo^j+ijoo T Oir ), 


(4.31) 


Therefore the matrix ^(p) remains bounded as r —>■ oo, that is, p —)■ 0. 

As p —)■ oo, we have r Vh- Using the boundary conditions fl2.25p . we hnd that 


dr,u:j = 0{r-rh), Q = S{rh)p\rh) + 0{r - th), (4.32) 


and that the entries of the matrix Cf are all 0(r — Vh) as r ^ Vh, so that £(p) —?■ 0 as p —)• oo, 
satisfying fl4.23p . This also means that the matrix £(p) is uniformly bounded on [0, oo). 

To check whether fl4.24l) is satished, we hrst note that, for r ~ r/j. 


p = -r* ~ -ph In (r - Vh ), 


(4.33) 
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where ph is a positive constant. Therefore, as p —)■ cx), 


£(p) ~ O (e Ph 


(4.34) 


Therefore —)■ 0 as p —)■ cxo. Therefore there exists a pi snch that for all p > pi, 

1 < < -i (4.36) 


Hence we have satished fl4.24p with b = 1/8. 

Therefore we have cast the gravitational sector pertnrbation eqnations in the form fl4.22p 
reqnired for the application of the nodal theorem, and all the conditions reqnired by the 
theorem are satished. We comment that this case was simpler to deal with than the sitnation 


in Ref. 


19l . becanse in that paper terms arising from non-spherically symmetric pertnrbations 


cannot be inclnded in Cl(p). 


Solitons 


For solitons, the system of gravitational pertnrbation eqnations needs fnrther transfor¬ 
mation before it is in the form reqnired for the application of the nodal theorem. 


Following Ref. Il9|, we dehne 


p = r 2 ^ 


(4.36) 


so that p —)■ 0 corresponds to r —)■ cx) and p —)• cx) corresponds to the origin. We also make 
a transformation of the pertnrbations: 


Soj = Xv, 


(4.37) 


where 

X = r2 {pS)~^ . 

The gravitational sector pertnrbation eqnations now take the form 

d?v 


where 


AX Xv = —— Mv, 

dp^ 


'—' . d f 1 dX\ 


(4.38) 


(4.39) 


(4.40) 
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Comparing with fl4.20p . as with the black hole case we take 


^(P) — 1'N-l- 


(4.41) 


To £x the matrices 5B(p) and £(p), we need to study the behaviour of the matrix M. as 

p —)■ 0 . 

As p —0, r —)■ cxD and 

\/3 _i ^ / _5\ \/3 


X = 


/o _1 ^ / _5\ VO 1 ^ 

-r + O [r 4 j = _^_p 2 + O ( p 2 1 . 

—A V / y—A V / 


\/—A 


Therefore, as p —>■ 0, 


d / 1 d:TA 3 


(4.42) 


(4.43) 


From the analysis of the previous subsection, we know that Ai = 0(1) as r ^ oo. This 
suggests that we should take 


»(rt = -p‘x4- () I„_„ €[p) = iX'M. 


dp dp 


With this choice, we have 


®( 0 ) = 


(4.44) 


(4.45) 


which is non-negative as required. 

To check the other conditions (14.231 I4.24p on the matrices ®(p) and ^(p), we need to 
examine their behaviour as p ^ oo, that is, r —)■ 0. In this case, using the boundary 
conditions fl2.23p 

1 *:^ / "7 \ 1 / •7 \ 

(4.46) 


1 3 f '^ \ 1 3 

X = —;=r4 -I- O ( r4 ) = —p=p~2 -p o (p~ 


V^o 


which gives 


3 


d / 1 dXA ^ -2 _i_ -4^ 


(4.47) 


Therefore 23(p) IXy-i as p —)■ cx) and the matrix f8(p) is uniformly bounded on [0, oo). 

We next turn to the behaviour of the matrix C(p) as p —?• 0 and r —)• oo. We have already 
seen that Ai = 0(1) as r —)■ oo. Then, using the dehnition of Ci(p) (I4.44p and the behaviour 
of X as p —)■ 0 (I4.42p . the matrix £(p) —)■ 0 as p —)■ 0. 

Therefore it remains to investigate the properties of (t(p) as p —)■ oo and r ^ 0. Using 
the boundary conditions fl2.23p . we first note that, as r ^ 0, 


dr^Oj = 0{r), Wj = 0{r^), 2 = ^ + 0 ( 1 ). 


(4.48) 


35 






















Therefore the behaviour of the entries of the matrix M fl4.9p as r —)■ 0 is: 

992 

+ (N-j-1) + 0 ( 1 ), 

Mj,k = 0(1). (4.49) 

Using the behaviour of X as p —)• cx) fl4.46p . we then have 

€(p) = 0(p~^) —)■ 0 as p —)■ cx), (4.50) 

so (I4.23P is satished. 

Looking at the remaining condition fl4.24l) . using the asymptotic forms fl4.46l IT~471 14.49p . 
we see that, as p —)■ 00 , 

iB(p) + p^€(p) = + 4£jv-i + 0(p-^) (4.51) 

where the matrix £]v-i is given by (14.141) . Since we know that the matrix £jv-i has only 
positive eigenvalues, we deduce that 53(p) + p^^(p) is positive for sufficiently large p and 
therefore fl4.24p is satisfied. 

Finally in this subsection we note that the eigenvalue problem we have in the gravitational 
sector fl4.39p is not exactly of the form (14.221) required for the application of the nodal 
theorem. However, since > 0 everywhere, this will not be a major difficulty in our 
analysis in section IIV El 

E. Existence of static solutions with no gravitational sector instabilities 

We are now in a position to prove, in this subsection, the existence of non-trivial (that is, 
non-embedded) su(iV) solitons and black holes which have no instabilities in the gravitational 
sector. For both solitons and black holes, our argument will use Amann and Quittner’s nodal 

theoreirt^Z. 

For black holes, in section IIV D 21 we have written the gravitational sector perturbation 
equations in the standard form (14.221) required for the application of the nodal theorem. In 
order to show that black hole solutions have no instabilities in the gravitational sector, it 
therefore suffices to show that the function ^(p) (14.271) has no zeros on an interval p G (c, 5), 
for small c and large h. 
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For soliton solutions, the argument is a little more involved. First of all, the gravitational 
sector perturbation equations fl4.39p take the form 

(fv 


0 Au = T)v = 


dp^ 


+ Mv, 


(4.52) 


where 0 = is a positive function, whereas the nodal theorem applies to the eigenvalue 
problem 'Dv = \v fl4.22p . Suppose that we are able to show that there exist su(A^) soliton 
solutions for which the function '^{p) fl4.27p has no zeros in the interval p G (c,f)). Then, 
applying the nodal theorem, the eigenvalue problem (14.221) has no negative eigenvalues. This 
means that the operator 2) is a positive operator. Then, if 2) is a positive operator, it must 
be the case that the eigenvalue problem (14.521) also cannot have any negative eigenvalues 
because 0 is a positive function. The upshot is that, for the soliton case as for the black hole 
case, if we can show that the function ^{p) (14.271) has no zeros in an appropriate interval, 
then there are no instabilities in the gravitational sector. 

I, we know that the equilibrium su(A^) solutions 


From the existence theorems in Ref. 


30 


of the held equations are analytic in r, A and the parameters at the origin or event horizon 
characterizing either soliton or black hole solutions. The matrices 21, 23 and 0 appearing in 
the operator 2) (I4.2np are analytic functions of the equilibrium held functions /i(r), S{r) and 
u)j(r) and r (and hence p) for values of p in our interval of interest (^). Standard existence 
theorems for ordinary diherential equations (see, for example. Ref. 1381) then tell us that the 
solutions Uj of the initial value problems (14.261) are also analytic functions of p, A and the 
initial parameters at either the origin or event horizon. Therefore, the function ^{p) (I4.27p 
is also analytic in p, A and the initial parameters at the origin or event horizon. 

In section UlI E 31 we proved the existence of embedded su(2) solitons and black holes for 
which Ci}(r)^ > 1 for all r. There we also showed that these embedded su(2) solutions have 
no instabilities in the sphaleronic sector. From section IIV C 31 they also have no instabilities 
in the gravitational sector. Pick such an embedded su(2) solution (either a soliton or a 
black hole). Fix c to be very small and c) to be very large. Then, from the nodal theorem, 
the function ^ will have no zeros on the interval (c, O) for this particular embedded su(2) 

, there exist genuinely su(A^) 


soliton or black hole. From the existence theorems in Ref. 


30 


solutions in a neighbourhood of this embedded 5u(2) solution. Since ^{p) is analytic in the 
parameters at the origin or event horizon which characterize the su(A^) solutions, providing 
the 5 u(A^) solutions are sufficiently close to the embedded 5u(2) solution, the function '^{p) 
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will continue to have no zeros in the interval (c,t)) for the su(iV) solitons or black holes. 

Therefore, if we consider su(A^) solutions sufficiently close to this stable embedded 5u(2) 
solution, using the nodal theorem (and considering the operator 'D for soliton solutions as 
described above), we have therefore proven that these su{N) solutions have no instabilities 
in the gravitational sector. 

In section IIII FI we showed that su(iV) solutions in a neighbourhood Tli of the above 
stable embedded su(2) solution have no instabilities in the sphaleronic sector. Having, in 
the current section, shown that su(iV) solutions in another neighbourhood 912 of the above 
embedded su(2) solution have no instabilities in the gravitational sector, we can deduce that 
those su(iV) solutions in the intersection of 9ti and 912 are stable under linear, spherically 
symmetric perturbations. 


V. CONCLUSIONS 


In this paper we have proven the existence of non-trivial, purely magnetic, spherically 
symmetric, su(A^) Einstein-Yang-Mills solitons and black holes in asymptotically anti-de Sit¬ 
ter space (with a negative cosmological constant A) which are stable under linear, spherically 
symmetric perturbations. 

The equilibrium solutions we consider are purely magnetic and spherically symmetric and 
the Yang-Mills held is described by Y — 1 functions 0Jj{r). With an appropriate choice of 
gauge, the perturbation equations for linear, spherically symmetric, perturbations decouple 
into two sectors: the sphaleronic sector and the gravitational sector. The sphaleronic sector, 
involving only gauge held perturbations, is easier to analyze and is considered in section 
m We hnd a series of inequalities fl3.33p on the equilibrium functions 0Jj{r) which are 
sufficient for there to be no instabilities in the sphaleronic sector. We hrst proved the 
existence of embedded su(2) solutions which satisfy these inequalities, before showing that 
su{N) solutions in a neighbourhood of these stable embedded su(2) solutions also have no 
instabilities in the sphaleronic sector. 

The gravitational sector is studied in section IIVI The metric perturbations can be elim¬ 
inated to leave a set of equations for gauge held perturbations. Our approach to proving 
stability in this sector follows Ref. [lOl, making use of a nodal theorem for a multidimensional 


Schrodinger system^. Again we can prove the existence of su{N) solutions, in a neighbour- 
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hood of stable embedded su(2) solutions, which have no instabilities in the gravitational 
sector. 

A natural question is how the stable su(A^) EYM black holes whose existence we have 
proven in this paper £t into the context of the “no-hair” conjecture as formulated by Bizort^: 


Within a given matter model, a stable stationary black hole is uniquely de¬ 
termined by global charges. 


It is argued in Ref. l40| that, for sufficiently large |A|, there exist N — 1 non-Abelian magnetic 


global charges which uniquely characterize 5u{N) EYM black holes, at least for large event 
horizon radius Vh and in a region of the parameter space which contains embedded su(2) 
black holes. For both the sphaleronic and gravitational sectors, our proof of the existence 
of stable su(iV) EYM black holes (and solitons) is valid for large |A|. Combining our results 
in this paper with those in Ref. 1^, we have evidence that at least some large stable su(Y) 
EYM black holes are uniquely determined by global charges, in accordance with Bizon’s 




“no-hair” conjecture (see also Ref. 

In this paper we have considered only purely magnetic, spherically symmetric solitons 
and black holes. The existence of su(iV) purely magnetic topological black holes has been 
proveii^, and solutions found numerically for the su(3) gauge group^. Very recently it has 
been shown that the argument we have presented here can be extended to show the stability 
of some of these su(iV) purely magnetic topological black holes^. Dyonic solitons and black 
holes in su(2) EYM have been found numerically^ and the existence of solutions where both 
the electric and magnetic gauge field functions have no zeros has been proveii^. Dyonic 
solutions have also been found numerically for the larger gauge group 5u(3}^. Recently the 
existence of dyonic soliton and black hole solutions of the su(Y) field equations has been 
proveii^. The existence of stable su(2) dyonic solutions has been proven very recently^ and 
it would be interesting to investigate whether our results in this paper on the existence of 
stable purely magnetic solitons and black holes in su{N) EYM in anti-de Sitter space can 
be extended to dyonic solutions. 

Finally, we comment that in this paper our focus has been the classical stability of su(Y) 
purely magnetic EYM black holes and solitons in anti-de Sitter space. We have considered 
only linear, spherically symmetric perturbations. The extension of our results to general 


linear perturbations is likely to be extremely challenging technically (see Refs. M and 
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for the 5u(2) case) and we would expect that at least some of the solutions which are stable 
under spherically symmetric linear perturbations will remain stable when general linear 
perturbations are considered. Going beyond classical stability, recent work has considered 
the thermodynamics of purely magnetic 5u(2) EYM black holes in anti-de Sitter spaced 


(see also Refs. 


22 


and 


401 ). In the su(2) EYM case, for generic (non-integer) magnetic 


charge there are two branches of asymptotically anti-de Sitter black holes, one of which is 


thermodynamically stable. It would be interesting to extend the work of Ref. |49 
larger su(iV) gauge group. 


to the 
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